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There  are  strong  needs  for  productive/quality  machining  strategies  of  notoriously  "difficult-to-machine"  aerospace 
materials,  automobile  parts,  surgical  instruments  and  super  alloys.  The  newest  EDM  technology  may  be  able  to  machine 
intricate  and  complex  3D  structures  and  also  circumvent  the  problems  encountered  in  mechanical  machining  methods. 
This  study  is  based  on  parametric  optimization  of  powder-mixed  EDM  machining  of  EN-31  tool  steel.  Copper  is  used  as 
micro-tool.  GRA  and  TOPSIS  are  used  for  optimization  of  parameters.  The  optimum  parametric  conditions  were  found  to 
be  150us/0.7DC/12A/2g/l.  Peak  current  (Ip)  was  found  to  be  the  most  influencing  factor  followed  by  cone  (D)  then 
time-on-pulse  (Ton)and  finally  duty  cycle  DC.  It  is  also  observed  that  with  powder-mixed  dielectrics.  MRR  was  found  to  be 
high  and  Ra  to  be  low.  Moreover,  the  machined  surface  develops  high  resistance  to  corrosion  and  abrasion. 

KEYWORDS:  Grey  Relational  Analysis,  Electro  Discharge  Machining 


Electro  Discharge  Machining  is  one  of  the  advanced  machining  process  where  the  material  removal  occurs  by 
melting  and  vaporization  due  to  the  thermal  energy  produced  between  the  work  piece  and  the  electrode.  All  electrically 
conductive  material  can  be  machined  by  EDM.  Powder  mixed  EDM  (PMEDM)  has  a  different  machining  mechanism  from 
the  conventional  EDM  [1].  In  this  process,  a  suitable  material  in  the  powder  form  is  mixed  into  the  dielectric  fluid  either  in 
the  same  tank  or  in  a  separate  tank.  For  better  circulation  of  the  powder  mixed  dielectric,  a  stirring  system  is  employed. 
For  constant  reuse  of  powder  in  the  dielectric  fluid,  a  modified  circulation  system  is  used.  The  experimental  setup  consists 
of  a  transparent  bath  like  container,  called  machining  tank.  It  is  placed  in  the  work  tank  of  EDM  and  the  machining  is 
performed  in  this  container.  To  hold  the  work  piece,  a  work  piece  fixture  assembly  is  placed  in  it.  The  machining  tank  is 
filled  up  with  dielectric  fluid.  To  avoid  particle  settling,  a  stirring  system  is  incorporated.  A  small  dielectric  circulation 
pump  is  used  for  proper  circulation  of  the  powder  mixed  dielectric  fluid  into  the  discharge  gap.  The  pump  and  the  stirrer 
assembly  are  placed  in  the  same  tank  in  which  machining  is  performed.  The  distance  between  powder  mixed  dielectric 
suction  point  and  nozzle  outlet  is  kept  as  short  as  possible  (10  in.)  in  order  to  ensure  the  complete  suspension  of  powder  in 
the  discharge  gap.  Magnetic  forces  were  used  to  separate  the  debris  from  the  dielectric  fluid.  For  this  purpose, 
two  permanent  magnets  areplaced  at  the  bottom  of  machining  tank.  The  various  powders  that  can  be  added  into  the 
dielectric  fluid  are  aluminum,  chromium,  graphite,  silicon,  copper  or  silicon  carbide,  etc.  The  spark  gap  is  filled  up  with 
powder  particles.  When  voltage  is  applied  the  powder  particles  get  energized  and  behave  in  a  zigzag  fashion  [2]. 
These  charged  particles  are  accelerated  by  the  electric  field  and  act  as  conductors. 
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Servo  control 


Figure  1:  PMEDM  Set-Up 

The  conductive  particles  promote  breakdown  in  the  gap  and  increase  the  spark  gap  between  tool  and  the 
workpiece.  Under  the  sparking  area,  the  particles  come  close  to  each  other  and  arrange  themselves  in  the  form  of  chain  like 
strutures  between  both  the  electrodes.  The  interlocking  between  the  different  powder  particles  occurs  in  the  direction  of 
flow  of  current.  The  chain  formation  helps  in  bridging  the  discharge  gap  between  both  the  electrodes.  Due  to  bridging 
effect,  the  insulating  strength  of  the  dielectric  fluid  decreases.  The  easy  short  circuit  takes  place,  which  causes  early 
explosion  in  the  gap.  As  a  result,  a  'series  discharge'  starts  under  the  electrode  area.  The  faster  sparking  within  a  discharge 
takes  place  causing  faster  erosion  from  the  workpiece  surface  and  hence  the  material  removal  rate  (MRR)  increases. 
At  the  same  time,  the  added  powder  modifies  the  plasma  channel.  The  plasma  channel  becomes  enlarged  and  widened  [3]. 
The  sparking  is  uniformly  distributed  among  the  powder  particles,  hence  electric  density  of  the  spark  decreases. 
Due  to  uniform  distribution  of  sparking  among  the  powder  particles,  shallow  craters  are  produced  on  the  workpiece 
surface.  This  results  in  improvement  in  surface  finish[4].In  this  process,  an  appropriate  kind  of  powder  is  mixed  into  the 
dielectric  fluid  when  a  voltage  of  80-320  V  is  applied  to  both  the  electrodes;  an  electric  field  in  the  range  of  105-107  V/m 
is  created.  The  spark  gap  is  filled  up  with  additive  particles  and  the  gap  distance  between  tool  and  workpiece  increased 
from  25-50  urn  to  50-150  urn  [5]. 

EXPERIMENTAL  DESIGN  AND  DATA  COLLECTION 

This  study  is  based  on  secondary  data  provided  from  the  research  done  by  H.K.  Kansal,  Sehijpal  Singh, 
P.  Kumar  [6].  They  optimized  the  process  parameters  of  powder  mixed  electrical  discharge  machining  (PMEDM). 
Response  surface  methodology  (RSM)  was  used  to  plan  and  analyze  the  experiments.  The  process  was  studied  with  a 
standard  RSM  design  called  central  composite  design  (CCD).  A  total  of  30  experiments  were  conducted  at  the  stipulated 
conditions.  In  this  design,  the  factorial  design  is  full  factorial  design  with  all  combinations  of  the  factors  at  the  two  levels 

Machining  Conditions 

The  machining  conditions  are  given  in  Table  1 . 
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Table  1:  Machining  Conditions 


Conditions 

Descriptions 

Workpiece  (+) 

EN-31  tool  steel 

Electrode  (-) 

Copper  microtool,  25mm  diameter 

Dielectric  fluid 

Kerosene. 

Powder 

Silicon  (Si) 

Powder  size  (um) 

20-30  um 

Peak  current  (A) 

3-12 

Pulse-on-time  (lis) 

50-150 

Duty  factor  (%) 

0.7-0.9 

Experimental  Matrix  and  Results  for  the  PEDM  Performance  Characteristics 


Table  2:  Design  of  Experimental  Matrix  [6] 


Expt  No 

Process  Parameters 

Pulse  on  Time,Ton  (us) 

Duty  Cycle  DC 

Peak  Current  ID  (A) 

Concentration  D  (g/1) 

1 

50 

0.7 

3 

0 

2 

150 

0.7 

3 

0 

3 

50 

0.9 

3 

0 

4 

150 

0.9 

3 

0 

5 

50 

0.7 

12 

o 

6 

150 

0.7 

12 

o 

7 

SO 

0  Q 

n 
yj 

150 

0.9 

12 

o 

9 

50 

0.7 

3 

2 

10 

150 

0.7 

3 

2 

11 

50 

0.9 

3 

2 

12 

150 

0.9 

3 

2 

13 

50 

0.7 

12 

2 

14 

150 

0.7 

12 

2 

15 

50 

0.9 

12 

2 

16 

150 

0.9 

12 

2 

17 

50 

0.8 

7.5 

1 

18 

150 

0.8 

7.5 

1 

19 

100 

0.7 

7.5 

1 

20 

100 

0.9 

7.5 

1 

21 

100 

0.8 

3 

1 

22 

100 

0.8 

12 

1 

23 

100 

0.8 

7.5 

0 

24 

100 

0.8 

7.5 

2 

25 

100 

0.8 

7.5 

1 

26 

100 

0.8 

7.5 

1 

27 

100 

0.8 

7.5 

1 

28 

100 

0.8 

7.5 

1 

29 

100 

0.8 

7.5 

1 

30 

100 

0.8 

7.5 

1 

The  ranges  of  these  parameters  were  selected  on  the  basis  of  preliminary  experiments  conducted  by  using  one 
variable  at  a  time  approach.  Table  3gives  the  levels  of  various  parameters  and  their  designation.  The  response  parameter  in 
the  present  study  was  MRR  and  surface  roughness  (SR).  In  each  test,  the  MRR  was  calculated  by  the  weight  loss  method. 
The  SR  was  measured  in  terms  of  arithmetic  mean  roughness  of  the  evaluated  roughness  profile  (Ra.  in  um)  by  using  a 
Surfcoder  SE1200  surface  testing  analyzer. 
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Table  3:  Design  Scheme  of  Process  Parameters  and  Their  Levels  [6] 


Factor  Symbol 

Parameter 

Levels 

Low  (-1) 

High  (+1) 

A 

Pulse  on  time,  Ton  (|is) 

50 

150 

B 

Duty  cycle,  DC 

0.7 

0.9 

C 

Peak  current,  Ip  (A) 

3 

12 

D 

Concentration,  Cone,  (g/1) 

0 

2 

Table  4:  Experimental  Matrix  with  Process  Parameters  and  Their  Levels 


1 
1 

i 
i 

i 
i 

i 
i 

o 

Z 

i 
i 

i 
i 



D 

i 
i 

i 
i 



A 
4 

D 

i 
i 



J 

1 
1 

1 
1 

0 

1 
1 

3 



7 
/ 

1 



o 
O 



Q 

1 

1 

1 

D 

1  n 

1 
1 

1 
1 

J 

1  1 
1 1 

1 
1 

1 
1 

J 

1Z 

1 

D 

1  ^ 

1 J 

1 

1 

14 

3 

1 

3 

3 

15 

1 

3 

3 

3 

16 

3 

3 

3 

3 

17 

1 

2 

2 

2 

18 

3 

2 

2 

2 

19 

2 

1 

2 

2 

20 

2 

3 

2 

2 

21 

2 

2 

1 

2 

22 

2 

2 

3 

2 

23 

2 

2 

2 

1 

24 

2 

2 

2 

3 

25 

2 

2 

2 

2 

26 

2 

2 

2 

2 

27 

2 

2 

2 

2 

28 

2 

2 

2 

2 

29 

2 

2 

2 

2 

30 

2 

2 

2 

2 

Table  5:  Observed  Values  of  Output  Responses  [6] 


MRR(mmVmin) 

Ra(um) 

I 

II 

1.84 

2.22 

2.03 

2.48 

1.96 

2.22 

1.98 

2.84 

2.41 

1.98 

2.30 

2.14 

2.19 

2.00 

2.10 

2.30 

2.33 

2.32 

2.88 

1.94 

2.41 

2.41 

1.88 

2.14 

22.30 

23.10 

22.70 

7.21 

6.42 

6.82 

23.41 

24.20 

23.81 

6.14 

9.36 

7.75 

21.65 

23.70 

22.68 

6.50 

8.91 

7.70 

Impact  Factor  (JCC):  1.6842 
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Table  5:  Contd. 


24.21 

25.20 

24.71 

6.36 

7.13 

6.75 

3.17 

3.17 

3.17 

2.11 

2.34 

2.22 

3.30 

3.47 

3.38 

1.21 

1.48 

1.34 

3.16 

3.24 

3.20 

4.68 

3.54 

4.11 

3.19 

3.56 

3.38 

1.31 

1.37 

1.34 

26.98 

25.21 

26.10 

5.5 

5.48 

5.50 

28.38 

29.44 

28.91 

6.19 

4.05 

5.12 

26.00 

26.39 

26.19 

6.40 

3.84 

5.12 

28.21 

29.20 

28.71 

5.00 

5.23 

5.12 

4.75 

4.63 

4.69 

3.55 

7.81 

5.68 

8.20 

6.19 

7.19 

7.39 

6.51 

6.95 

7.19 

5.20 

6.21 

6.00 

6.02 

6.01 

7.43 

7.14 

7.29 

7.32 

5.47 

6.40 

2.10 

1.95 

2.03 

1.71 

1.89 

1.80 

25.50 

24.52 

25.01 

6.91 

5.06 

5.99 

5.20 

7.19 

6.21 

7.55 

6.07 

6.81 

6.94 

7.65 

7.29 

4.20 

5.36 

4.78 

6.18 

6.21 

6.19 

7.32 

5.47 

6.39 

6.70 

7.28 

6.99 

6.30 

7.30 

6.80 

6.22 

8.09 

7.15 

6.51 

7.09 

6.80 

7.93 

6.59 

7.26 

6.30 

6.09 

6.19 

6.61 

7.15 

6.88 

6.75 

5.55 

6.15 

7.58 

6.82 

7.20 

5.11 

6.85 

5.98 

Grey  Relational  Analysis 

The  grey  relational  analysis  (GRA)  is  one  of  the  powerful  and  effective  soft-tool  to  analyze  various  processes 
having  multiple  performance  characteristics.  Generally,  GRA  is  carried  out  for  solving  complicated  problems  which  have 
interrelationships  among  the  designated  performance  characteristics [7]. 

Normalization  of  the  Experimental  Results  in  GRA 

Normalization  is  the  process  in  which  transformation  of  input  data  takes  place  to  an  evenly  distributed  data  in  a 
scale  range  between  0  and  1.  The  experimental  results  for  the  responses  i.e.  MRR  and  Ra  are  normalized  using  equations 
(l)-(2)  where  Xy  is  the  normalized  value  of  y^f  or  response  j  (j  =  1,2,3,  of  experiment  I  (/  =  1,2,3,. .  .,m)  [8]. 

If  the  response  is  of  higher-the-better  type,  then 

x(/=(^-min(^))/(max(^)-min(^))  (D 

If  the  response  is  of  smaller-the-better  type,  then  normalized  value  jc#  is  expressed  as 

xij=(max(yij)-yij)/(max(yij)-win(yij))  (2) 

If  the  normalized  value  xtj  for  a  response  j  of  experiment  i  is  equal  to  1  or  nearer  to  1,  then  it  is  said  that  the 
performance  of  that  particular  experiment  i  is  best  for  the  response  j.  That  normalized  value  is  termed  as  reference  value 
(x0j)  for  7'th  response  [9] . 

In  the  present  research  study,  MRR,  is  of  higher-the-better  type  response  whereas  Ra  is  smaller-the-better  type 
response. 
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Table  6:  Normalization  Values  of  Results 


1 

0 

0.862715 

2 

0.014136905 

0.875195 

3 

0.002604167 

0.847114 

4 

0.014136905 

0.875195 

5 

0.768973214 

0.145086 

6 

0.810267857 

0 

7 

0.768229167 

0.0078 

8 

0.84375 

0.156006 

9 

0.042410714 

0.862715 

10 

0.050223214 

1 

11 

0.043526786 

0.567863 

12 

0.050223214 

1 

13 

0.89546131 

0.351014 

14 

1 

0.410296 

15 

0.898809524 

0.410296 

16 

0.992559524 

0.410296 

17 

0.098958333 

0.322933 

18 

0.191964286 

0.124805 

19 

0.155505952 

0.271451 

20 

0.195684524 

0.210608 

21 

0 

0.928237 

22 

0.854910714 

0.274571 

23 

0.155505952 

0.146646 

24 

0.195684524 

0.463339 

25 

0.154761905 

0.212168 

26 

0.18452381 

0.148206 

27 

0.19047619 

0.148206 

28 

0.194568452 

0.24337 

29 

0.180431548 

0.24961 

30 

0.19233631 

0.276131 

Table  7:  Deviation  from  Ideal  Sequence  of  Ideal  Solution 


Deviation  of  MRR 


Deviation  of  Ra 


1 

1 

0.137285 

2 

0.985863 

0.124805 

3 

0.997396 

0.152886 

4 

0.985863 

0.124805 

5 

0.231027 

0.854914 

6 

0.189732 

1 

7 

0.231771 

0.9922 

8 

0.15625 

0.843994 

9 

0.957589 

0.137285 

10 

0.949777 

0 

11 

0.956473 

0.432137 

12 

0.949777 

0 

13 

0.104539 

0.648986 

14 

0 

0.589704 

15 

0.10119 

0.589704 

16 

0.00744 

0.589704 

18 

0.808036 

0.875195 

19 

0.844494 

0.728549 

20 

0.804315 

0.789392 

Impact  Factor  (JCC):  1.6842 
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Table  7:  Contd. 


21 

1 

0.071763 

22 

0.145089 

0.725429 

23 

0.844494 

0.853354 

24 

0.804315 

0.536661 

25 

0.845238 

0.787832 

26 

0.815476 

0.851794 

27 

0.809524 

0.851794 

28 

0.805432 

0.75663 

29 

0.819568 

0.75039 

30 

0.807664 

0.723869 

Grey  Relational  Coefficient 

The  grey  relational  coefficient  is  calculated  to  determine  the  closeness  of  Xy  to  x0j.  Higher  value  of  grey  relational 
coefficient  means  Xy  is  closer  to  x0j.  Grey  relational  coefficient  is  calculated  based  on  equation  (3)  [8]. 

?(^^=(A^   ,mandj  =  1,2,3,.  ...n  (3) 

where,  A;  =\X0j  ~Xij\ 

Arin  =irin{^,i  =lZ%.-.,n%  7=1,2,3,....,^} 

□  =  distinguishing  coefficient,  iD(0,l) 

The  normalized  values  of  each  of  the  responses  for  all  30  experiments  are  used  to  calculate  the  grey  relational 
coefficient  using  equation  (3).  The  distinguishing  coefficient  □  is  taken  as  0.5. 

Grey  Relational  Grade 

Grey  relational  grade  is  the  weighted  sum  of  the  grey  relational  coefficients  for  a  particular  experiment  and  it  is 
calculated  using  equation  (4)  [9] . 

n 

r(XQ,X.  )  =  Y,WJ'  y(x0j>xij  )  for  1  =  h2,3,...,m  and;  =  1,2,3,....,^  (4) 
where,  D(Xo,  xi)  is  grey  relational  grade  between  comparability  sequence  Xt  and  reference  sequence  X0  and 
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Table  8:  Grey  Relational  Grades 


Expt  No. 

MRR(mg/min) 

Ra(um) 

Coeff  of 
MKK(mg/min ) 

Coeff  of  Ra(um) 

Grey  Relational 
Grades 

Order 

1 

1 

0 137285 

vj  •  x  *y  /  jl*  vj  +j 

0  333333 

vy  •  *y  *y  *y  *y  *y  *y 

0  784578 

vy  •  /  w  r_/  /  vj 

0  55896 

vj  •        yj  y  vj 

11 

2 

0  985863 

V7  .  y  U  ^7      V7  ^7 

0 1 24805 

0  336505 

Vy.^7^7V7^7Vy^7 

0  80025 

v7 .  u  vy  vyz-^y 

0  56838 

vy .  ^7  vy  vj  ^7  vj 

8 

3 

0  997396 

vy  •  y  y  i  *y  y  vy 

0 152886 

VJ  •  X  +J  A^yJ  yj  VJ 

0  333913 

vy.^y^y^yv  i 

0  76583 

vj  •  i  vj  +j  kj  *y 

0.54987 

12 

4 

0  985863 

0 124805 

0  336505 

vy  •  *y  *y  vj^j  vj \j 

0  80025 

VJ  •  VJ  VJ vj  ^  +j 

0  56838 

VJ  •  +J  VJ  yj  *y  yj 

8 

5 

0  231027 

0  854914 

V7  .  L7  ^7    1  y  1  1 

0  683969 

V7  .  V7  L7  ^7  y  V7  V 

0  369027 

V7  .  ^'  V7  V  V7  / 

0.5265 

13 

6 

0 189732 

V7 .  1  U  y  /  ~J 

1 

0.724919 

0  333333 

vy.^y^y^y^y^y^y 

0.52913 

11 

7 

0.231771 

0.9922 

0  683274 

0  335076 

vj  •  *y  *y  +j  vj  i  vj 

0  50918 

vj  •  +j  vj  y  l  yj 

12 

8 

0 15625 

V7  .  1  ^7  UAJ 

0  843994 

\7.U~J77~ 

0  761905 

V7  .  /  V7  1  y  V7  ^7 

0  372025 

0  56697 

vy .  ^7  vy  vy  v  / 

9 

9 

0  957589 

V7  .  V  ^7  /  ^7  U  7 

0 137285 

V7  .  1  ^7  /       VJ  -7 

0  343032 

vy.^7i^7vy^7z- 

0  784578 

W.  /  VJ~-7  /  VJ 

0  56381 

vy .  ^7  vy  ^y  vj  i 

10 

10 

0.949777 

o 

0  344881 

vy  •  *y  i  r  vj  u  i 

1 

0.67244 

3 

11 

0  956473 

vy  •  y    u  i  /  w-/ 

0  432137 

vj  •  r    ^  i  / 

0  343295 

V7  .       1          y  ^7 

0  536402 

vj  •    *y  vj  vvj  ^ 

0  43985 

V7  .   1  ^7  y  (_  )  ^7 

14 

12 

0.949777 

o 

0  344881 

V7 .         ru  vj  i 

1 

0.67244 

3 

13 

0 104539 

V7  .  1  V/  I  J  ^7  V 

0  648986 

V7 .  V7  l  Uy  UU 

0  827076 

vy.vjz-/vy/vy 

0  435166 

V7  .   1  ^7  ^7  1  V7  V7 

0.63112 

5 

14 

o 

0  589704 

vj  •    kj  y  i  vj r 

1 

0  45884 

vj  •  r    vj  vj  r 

0.72942 

1 

15 

0.10119 

0  589704 

Ut^uy  /  V7  i 

0  831684 

vy .  U  ^7  1  V7  L7  1 

0  45884 

vy  V7V7~ 

0.63112 

4 

16 

0.00744 

0  589704 

Uijuy  /  V7  i 

0  985338 

Vy.VVJ_7.-7.-7VJ 

0  45884 

W  .~_7  V7V7~ 

0  72209 

vy .  /  Z/A/Uy 

2 

17 

0  901042 

vy  •  y  vj  i  vy 

0  677067 

W  •  V/  /    /  WVy  / 

0  356877 

vy  •  *y     VJ  VJ  1  1 

0  424785 

vj •         r  / 

0  39083 

vy  •  ■^j  y  vj  yj  *y 

17 

18 

0  808036 

W  •  VJ  \J\j\J *J  VJ 

0  875195 

vj  •  yJ  i  +j  x  y  +j 

0  382252 

v./  •  *y  vj    ;L*  *y  jL* 

0  363585 

w  •  *y  vj  *y    vj  *y 

0.37292 

25 

19 

0.844494 

0  728549 

V7  .  /       K )  ^7   1  y 

0  371887 

0  406984 

V7  .   1  V7  V7  y  V)  1 

0.38944 

19 

20 

0  804315 

vj  •  yj  vj r     l  +j 

0  789392 

vy  •  i  vj  y  *y  y  ^ 

0  383343 

vj  •  *y  vj  *y  *y  r  *y 

0  38778 

vy  •  *y  vj  1  i  yj 

0  38556 

vy  •  ju  j  j  vy 

21 

21 

1 

0  071763 

vy  •  vy  /±/  vj *y 

0  333333 

vy  •  *y  *y  *y  *y  *y  *y 

0  874488 

vj  •  yj  I    i  r  yj  yj 

0  60391 

vj  •  vy  vy  *y  y  i 

6 

22 

0 145089 

0  725429 

V7  .  /  Z  JTZj7 

0  775087 

vy .  /  /  — 7  vy  VJ  / 

0  40802 

vy  .~vy  vjvyzw 

0  59155 

vy .  ^7  V  J-  ^7  ^7 

7 

23 

0.844494 

0  853354 

0  371887 

0  369452 

V7  .  ^7  vs  y  1  ^7 

0  37067 

vy .  ^y  /  vy  vy  / 

26 

24 

0.804315 

0.536661 

0.383343 

0.482318 

0.43283 

15 

25 

0.845238 

0.787832 

0.371681 

0.388249 

0.37997 

22 

26 

0.815476 

0.851794 

0.380091 

0.369879 

0.37499 

24 

27 

0.809524 

0.851794 

0.381818 

0.369879 

0.37585 

23 

28 

0.805432 

0.75663 

0.383015 

0.39789 

0.39045 

18 

29 

0.819568 

0.75039 

0.378912 

0.399875 

0.38939 

20 

30 

0.807664 

0.723869 

0.382361 

0.40854 

0.39545 

16 

Table  9:  Average  and  Range  of  Grey  Relational  Grade 


Average  Grey  Relational  Grade 

Max-Min 

Pulse  on  time,Ton  (urn) 

0.533471 

0.423338 

0.600239 

0.176901 

Duty  cycle  DC 

0.574356 

0.422401 

0.560607 

0.151955 

Peak  current  Ip  (A) 

0.57756 

0.38736 

0.60412 

0.21676 

Concentration  D  (g/1) 

0.52756 

0.420026 

0.610569 

0.190543 

Mean  Value  of  the  Grey  Relational  Grade  =  0.18404. 

Impact  Factor  (JCC):  1.6842 


Index  Copernicus  Value  (ICV):  3.0 
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Experiment  number 

Figure  2:  Plot  of  Experiment  Number  vs  Grey  Relational  Grades 

The  above  graph  is  plotted  with  software  Origin  The  experiment  which  has  highest  grey  relational  grade  is  said  to 
be  best  choice  of  all  the  runs.  Table  12  shows  that  experiment  number  14  is  the  best  run.  The  optimized  process  parameters 
are  150us/0.7DC/12A/2g/l.  Ip  was  found  to  be  the  most  influencing  factor  followed  by  cone  D  and  then  Ton?  and  finally  and 
the  least  influencing  factor  DC. 

TOPSIS  Method 

TOPSIS  (Technique  for  Order  Preference  by  Similarity  to  Ideal  Solution)  is  one  of  the  many  multi  criteria 
decision  making  method(MCDM)  and  was  firstly  proposed  by  Hwang  and  Yoon  [10].  The  basic  concept  of  this  method  is 
that  the  chosen  alternative  (appropriate  alternative)  should  have  the  shortest  distance  from  the  positive  ideal  solution  and 
the  farthest  distance  from  negative  ideal  solution.  Positive  ideal  solution  is  a  solution  that  maximizes  the  benefit  criteria 
and  minimizes  adverse  criteria,  whereas  the  negative  ideal  solution  minimizes  the  benefit  criteria  and  maximizes  the 
adverse  or  cost  criteria. 


A 


Attribute  Xi  (increasing  preference) 

Figure  3:  Topsis  Graph  [10] 

The  experimental  data  used  in  present  study  is  based  on  the  research  done  by  H.K.  Kansal,  Sehijpal  Singh, 
P.  Kumar  [6]. 

STEPS  FOR  TOPSIS 
Step  1:  Decision  Matrix 

The  first  step  is  to  formulate  decision  matrix  with 'm'  alternatives  and  'n'  attributes  [11].  In  this  case  study  there 
are  30  number  of  experiments  (alternatives)  and  two  output  responses  MRR  and  Ra  (attributes).  Hence  the  decision  matrix 
is  of  the  order  30  X  2. 


www.tjprc.org 
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4.1 


Am 


-X, 


!^ml     \m-  Xmi 


Figure  4:  Decision  Matrix 


Table  10:  Decision  Matrix  with  30  Alternative  and  3  Attributes 


MRJttniiu^/min)^ 

1 

2.03 

2.22 

2 

241 

214 

3 

2.10 

2.32 

4 

2.41 

2.14 

5 

22.70 

6.32 

6 

23.31 

7.75 

7 

22.63 

"."C 

S 

24.71 

6.75 

9 

3.17 

2.22 

10 

3.33 

134 

11 

3.20 

411 

12 

3.33 

1.34 

13 

2610 

5.50 

14 

23.91 

512 

15 

2619 

512 

16 

23.71 

5  12 

17 

469 

5.63 

IS 

7  19 

695 

19 

621 

601 

20 

7.29 

6.40 

21 

2.03 

1  SO 

22 

2501 

5.99 

23 

6.21 

6.31 

24 

7.29 

4.73 

25 

619 

6.39 

26 

6.99 

6.30 

27 

7.15 

6.SC 

23 

7.26 

619 

29 

6.33 

615 

30 

7.20 

5.93 

Step  2:  Normalization  Matrix 

Normalization  matrix  is  formed  in  order  to  transform  the  various  attribute  dimensions  into  non-dimensional 
attributes,  which  allows  comparison  across  the  attributes  [11-13].  It  is  found  out  by  using  the  following  equation. 


■(5) 


\4 


Impact  Factor  (JCC):  1.6842 


Index  Copernicus  Value  (ICV):  3.0 
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Table  11:  Normalized  Matrix 


MRRtmm-.-'min)  \;i 

i 

2.03 

2.22 

2 

2.41 

2.14 

3 

2.10 

2.32 

4 

2.41 

2.14 

5 

22.70 

6.32 

6 

23. SI 

7.75 

7 

22.63 

7.70 

S 

24.71 

6.75 

9 

3.17 

2.22 

10 

3  J8 

1.34 

11 

3.20 

4.11 

12 

3.33 

1.34 

13 

26.10 

5.50 

14 

23.91 

5.12 

15 

2619 

5  12 

16 

23.71 

5  12 

17 

469 

5.63 

IS 

7.19 

695 

19 

621 

601 

20 

7.29 

640 

21 

203 

1  so 

22 

25  01 

5.99 

23 

621 

681 

24 

7.29 

478 

25 

619 

639 

26 

699 

630 

27 

7.15 

6  SO 

23 

126 

6  19 

29 

6.33 

6.15 

30 

7.20 

5.9S 

Step  3:  Weighted  Normalized  Decision-Making  Matrix 

Multiply  each  column  (ry)  by  weights  'wj'  assigned  to  attributes  to  get  Vy  [84-87]. 

vu-  =  wj*rij 

In  this  case  equal  weights  are  given  to  both  attributes  MRR  and  Ra 


Table  12:  Weighted  Normalized  Matrix 


iJ 

0.5 

tat  No 

Ra(Lim)  &h 

I 

0.012615 

2 

0.01497* 

0.03*144 

3 

0.013049 

0.0391 84 

4 

0.01497* 

0.03*144 

5 

0  14105* 

0.115187 

6 

0.147954 

0.130894 

7 

0.140932 

0.13005 

8 

0.15354* 

0.114005 

9 

0.01969B 

0.037495 

10 

0.021003 

0.022*32 

11 

0.019BB5 

0.0*9416 

12 

0.021003 

0.022*32 

13 

0.162184 

0.092893 

14 

0.179*45 

0.08*475 

15 

0.1*2743 

0.08*475 

16 

0.178402 

0.08*475 

17 

0.029144 

0.095933 

IS 

0. 044*7 B 

0.117383 

19 

0.03  85  89 

0  10150* 

20 

0.0453 

0.10B093 

21 

0.012*15 

0.030401 

22 

0.15541 

0  1011*9 

23 

0.038589 

0.115018 

24 

0.0453 

0.080732 

25 

0.038464 

0.107924 

2* 

0.04343* 

0.114B49 

27 

0.04443 

0.114849 

28 

0.045113 

0.104547 

29 

0.042752 

0.103871 

30 

0.04474 

0  101 
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Step  4:  Determine  Ideal  and  Negative  Ideal  Solution 

Ideal  solution; 

A+  =  {(max  Vijj  G  J),  (mfn  Vy j  £  J)  i  =  1,  2|  3,  m}[84-87]  (6) 

Negative  ideal  solution; 

A"  =  {(min  Vyj  G  J),  (ma|x  Vy j  G  J)i=  1,  2,|3,  m}  (7) 

where; 

J  =  {j  =  1,  2,  3,  n;  j  associated  with  benefit  criteria} 

J  =  {j  =  1,  2,  3,  n;  j  associated  with  cost  criteria} 

In  this  case; 
Benefit  criteria:  MRR. 
Cost  criteria:  Ra. 

From  weighted  normalized  decision-making  matrix,  ideal  solution  A*, 
A+  =  (0.179645;  0.022632) 

From  weighted  normalized  decision-making  matrix,  negative  ideal  solution  A", 
A-=  (0.012615;  0.130894). 
Step  5:  Determine  Separation  from  Ideal  Solution 

The  separation  measure  Si*  i.e.  the  distance  of  each  vy  from  ideal  solution; 
A+  =(0.179645;  0.022632).  is  calculated  by  using  the  following  equation  [12]. 
Si*  =  [  L  (vj*-  vy)2  ]  1/2for  each  row  j  (8) 

Table  13:  Separation  Measure  from  Ideal  Solution  Matrix 


Si* 

1 

2 

0. 165223 

3 

0167416 

0. 165223 

5 

0. 1002.7  7 

6 

0.112S05 

1 

O.H4131 

3 

0.09-5027 

S» 

0.160636 

io 

HI*  1  5  3-6-42 

11 

O  16647 

12 

0. 15 3642 

13 

0.0723-9-S 

14 

0.063S43 

15 

0.066042 

16 

C.Cl'S?:"5 

17 

0.167403 

13 

0.16+905 

1-P 

O 161611 

20 

0.159224 

2-1 

0_1672.11 

22 

0.032191 

23 

0.16S61S 

24 

0.1463  71 

25 

0„16494fi 

2-15 

0.16449 

27 

0.16366S 

2S 

0.15  7  503 

2-P 

O 1591S4 

30 

0.156015 

Impact  Factor  (JCC):  1.6842 


Index  Copernicus  Value  (ICV):  3.0 
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Step  6:  Determine  Separation  From  Negative  Ideal  Solution 

The  separation  measure  Sii.e.  the  distance  of  each  vy  from  negative  ideal  solution;A-  =  (0.012615;  0.130894),  is 
calculated  by  using  the  following  equation  [11,12]. 

Si-  =  [E(Vj-Vij)2f(9) 

Table  14:  Separation  Measure  from  Negative  Ideal  Solution  Matrix 


Sr 

1 

o.og-33-gi^ 

0.0947S 

3 

0 .09-1711 

4 

0.0947S 

5 

0.129-3-9'S 

<S 

0. 13  5339- 

T 

0.12S31P 

S 

0.14193P 

0.093<S<S7 

0.10S5S6 

11 

O -CW5190H5 

12 

0.10S5S6 

13 

0. 154321 

14 

0.172S35 

15 

0.15<S5-G1 

1-6 

O 171634 

17 

0.03S672 

IS 

O.O347-04 

19 

0.03P221 

2tJ 

0.03PS52 

21 

0.1O04-93 

22 

0.145S5G 

23 

0030441 

24 

0.059-371 

25 

0.0345S 

2G 

0.034747 

27 

0.03  5632 

25 

0.041S37 

29- 

0.04047S 

30 

0.043SS3 

Step  7:  Relative  Closeness  Coefficient 

The  relative  closeness/closeness  coefficient  (Q*)  to  the  ideal  solution  is  then  calculated  with  the  help  of  following 
equation  [12-13]. 

Q*  =  Si-/  (Si*  +Si-);  0<  Ci<l;  i=  1,  2,  3  m. 

Q*=  1;  if  Ai  =  A+ 

Q*  =  0;  if  Ai  =  A- 

Table  15:  Relative  Closeness  Coefficient  Matrix 


c  * 

0.35772S 

2 

0.364534 

3 

0.3  53523 

4 

0364534 

5 

0.51533M 

6= 

O  545404 

T 

0. 5251  52 

S 

0.59-se-se. 

5 

0.3<SS325 

IO 

0406343 

11 

0.271072 

12. 

0406343 

13 

0.-SSOI571 

14 

0.73025  5 

15 

0.703315 

H5 

0.72SS42 

17 

O. IS  76-55 

IS 

0. 1 7423 

15 

0. 155252 

20 

0  2001S4 

2.1 

0.3  7  53  SS 

22 

0.-6355SS 

23 

0. 152525 

24 

0.250255 

2J5 

0.173313 

2tf 

O  1744 

27 

0.17S7S4 

23 

0.205S71 

25 

0.20273 5 

30 

0.215526= 
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Step  8:  Average  Closeness  Coefficient 

Average  closeness  coefficient  is  the  weighted  sum  of  the  closeness  coefficients  for  a  particular  experiment  and  it 
is  calculated  using  equation  (5)  [11-13]. 

C°(xo,  x0=  -Y(x0j,  Xij)  for  /  =  l,2,3,...,m  and/' =  1,2,3,.  (10) 

where,  C(x0,  x{)  is  average  closeness  coefficient  between  comparability  sequence  Xt  and  reference  sequence  x0  and 

The  experiment  which  has  highest  relative  closeness  is  said  to  be  best  choice  of  all  the  runs.  The  values  of  average 
closeness  coefficients  and  range  are  shown  in  Table  16. 


Table  16:  Average  and  Range  of  Relative  Closeness  Coefficients 


Parameters 

Average  Relative  Closeness  Coefficients 

Max-Min 
(Range) 

Level  1 

Level  2 

Level  3 

0.228916139 

Pulse  on 
time,Ton  (um) 

0.446138778 

0.251025083 

0.479941222 

Duty  cycle  DC 

0.467994667 

0.467994667 

0.461817222 

0.006177445 

Peak  current  Ip 
(A) 

0.363132667 

0.196601167 

0.635512556 

0.438911389 

Concentration  D 
(g/D 

0.425620222 

0.2442475 

0.509496556 

0.265249056 

Mean  Value  of  the  relative  closeness  coefficients  =  0.288211271. 

a_  1  _|  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,— 

a  s  10  io  zi-  3<a 

Experiment  Num 


Figure  5:  Plot  of  Experiment  Number  vs  Relative  Closeness  Coefficient 

The  above  graph  is  plotted  with  software  Origin. 

By  using  Topsis,  it  was  found  that  the  experiment  number  14  has  the  maximum  value  of  grey  relational 
coefficient,  which  is  the  best  run.  The  optimized  process  parameters  i.e.  pulse-on-time  Ton,  duty  cycle  DC,  peak  current  \ 
and  concentration  D  were  found  to  be  150us/0.7DC/12A/2g/l  respectively.  Peak  current  Ip  was  found  to  be  the  most 
influencing  parameter  followed  by  concentration  D,  and  then  pulse-on-time  Ton5  and  finally  and  the  least  influencing  factor 
duty  cycle  DC. 


Impact  Factor  (JCC):  1.6842 


Index  Copernicus  Value  (ICV):  3.0 
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(c)  Effect  of  Ip  and  D  on  MRR  (f)  Effect  of  Ip  and  D  on  Ra 

Figure  6:  Surface  Plots  of  Various  Process  Parameters  vs  Output  Responses 

RESULTS  AND  DISCUSSIONS 

By  using  GRA,  it  was  found  that  the  experiment  number  14  has  the  maximum  value  of  grey  relational  coefficient, 
which  is  the  best  run  as  shown  in  Table  8.  The  optimized  process  parameters  were  found  to  be  150us/0.7DC/12A/2g/l 
respectively.  Ip,was  found  to  be  the  most  influencing  factor  since  it  has  the  maximum  value  of  range  of  average  grey 
relational  grades  as  shown  in  Table  15,  followed  by  cone  D  and  then  Tori5  and  finally  the  least  influencing  factor  DC. 

Ip>D>Ton>DC 
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With  Topsis,  experiment  number  14  was  found  to  be  the  optimum  run.  The  optimized  parameters  were  found  to 
be  150us/0.7DC/12A/2g/l.  Ip  was  found  to  be  the  most  influencing  parameter  followed  by  cone  D,  and  then  Ton^  and  finally 
and  the  least  influencing  factor  duty  cycle  DC. 

Ip>D>Ton>DC 

Analysis  of  Material  Removal  Rate 

It  was  observed  that  MRR  increased  with  increase  of  pulse-on-time  and  peak  current.  Surface  plot  indicates  that 
maximum  MRR  is  achieved  at  12  A  and  150  (is.  At  low  value  of  pulse-on-time  Ton,  high  MRR  is  achieved  at  high  value  of 
peak  current  Ip  (12A).  When  silicon  powder  was  suspended  into  the  dielectric  fluid  of  EDM,  an  enhanced  rate  of  material 
removal  rate  was  achieved.  This  is  due  to  the  fact  that  suspended  powder  causes  bridging  effect  between  both  the 
electrodes,  facilitates  the  dispersion  of  discharge  into  several  increments  and  hence  increases  the  MRR.  Maximum  MRR  is 
achieved  with  2g/l  concentration  of  powder.  Hence  more  improvement  in  MRR  is  expected  at  still  higher  concentration 
level  of  silicon  powder.  At  low  concentration  higher  MRR  is  achieved  at  high  peak  current  Ip  i.e.  12  A. 

Analysis  of  Surface  Roughness  Ra 

It  was  observed  that  surface  roughness,  Ra  increases  with  increase  of  pulse-on-time  Ton  and  peak  current  Ip.  This  is 
due  to  their  dominant  control  over  the  input  energy.  Surface  plot  showsminimumRa  is  achieved  at  Ton  150us  and  Ip  3  A  At 
high  values  of  Ton  minimum  Ra  is  achieved  at  low  value  peak  current.  Addition  of  silicon  powder  into  the  dielectric  fluid 
improved  surface  finish  of  the  workpiece.  Minimum  surface  roughness  is  achieved  at  high  concentration  2g/l.  Hence 
surface  roughness  is  expected  to  be  less  with  increasing  concentration  of  silicon  powder. 

CONCLUSIONS 

This  study  yielded  the  following  conclusions. 

•  MRR  and  surface  roughness  Ra  increases  with  increase  of  pulse-on-time  Ton,  and  peak  current  Ip 

•  An  enhanced  rate  of  material  removal  was  achieved   when  powder  abrasives  i.e.  silicon  powder  was  suspended 
into  the  dielectric  fluid. 

•  Surface  roughness  Ra  improves  with  addition  of  silicon  powder  and  is  expected  to  improve  more  with  more 
concentration  of  powder. 

•  GRA  and  Topsis  yielded  the  same  optimum  experiment  run  i.e.  14.  The  optimized  process  parameters 
are.l50us/0.7DC/12A/2g/l. 

•  With  both  GRA  and  Topsis,  peak  current  Ip  was  found  to  be  the  most  influencing  parameter  followed  by 
concentration  D,  and  then  pulse-on-time  Ton,  and  finally  and  the  least  influencing  factor  duty  cycle  DC. 
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